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I. INTRODUCTION 

The interplay between quantum effects and disorder is 
a subject actively studied for many decades, both theo- 
retically and experimentally. It plays a particularly im- 
portant role in condensed matter physics, where, in a first 
approximation, a crystal is modeled as independent elec- 
trons interacting with a perfectly periodic lattice. The pi- 
oneering works of Bloch and Zener [3, Q] showed however 
that most predictions based on this model are not verified 
in real crystals. For example, the Bloch theory predicts 
fully delocalized wavefunctions implying a ballistic trans- 
port of the electrons through the crystal. Moreover, in 
the presence of a constant bias potential, Zener predicted 
an oscillatory motion (the Bloch-Zener oscillations) due 
to quantum interference effects. This contradicts well- 
known experimental facts at least in usual conditions. 

An obvious possible explanation of these contradic- 
tions is the fact that there are no perfect crystals: In 
a real crystal some sites may be randomly occupied by 
ions of a different nature, thus breaking the periodicity of 
the lattice. In 1958, Anderson considered this approach 
and postulated that the dominant effect of the disorder 
is to change randomly the on-site energy. Starting from 
this assumption, he constructed a simple model [3| of 
a single-electron interacting with a lattice in the tight- 
binding approximation: 

#tb = ^2e jn \jn)(jn\ + ^ ^WmU^XM ■ (1) 

Here Cj n are the energies associated with the states la- 
beled by n at the sites j of the lattice, and the non- 
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diagonal elements Vj n ^^ denote the matrix elements be- 
tween these states. The diagonal part of the Hamiltonian 
corresponds to the potential energy and the non-diagonal 
part to the kinetic energy in a continuous space descrip- 
tion. Disorder is introduced by giving the site energies 
tj n a random distribution. Anderson thus showed that 
the electron wavefunctions can be localized by the disor- 
der. This is naturally in sharp contrast with the predic- 
tion of the Bloch model. 

The phenomenon of localization has its most striking 
manifestation in the transport properties of random me- 
dia. If particle-particle interactions are negligible, expo- 
nentially localized states cannot contribute to transport 
at zero temperature since the coupling to phonons is neg- 
ligible. Anderson localization as a consequence of the 
presence of disorder is one of the fundamental ingredi- 
ents for the understanding of the existence of insulators 
and metals, and, in particular, the transition between the 
insulating and the metallic states of matter. An insula- 
tor is associated with localized states of the system while 
a metal generally displays diffusive transport associated 
with delocalized states. 

It was later shown that the 3D Anderson model dis- 
played a phase transition between a localized and a 
diffusive phase, the so-called Anderson metal-insulator 
transition: If the disorder is below a critical level, 
the localization disappears and one recovers a metal- 
lic (conductor) behavior Q. The link between the 
disorder-induced metal-insulator transition and second- 
order phase-transitions was established by reformulating 
the problem in terms of the renormalization group 0, 
Based on Wegner's work and the ideas of Thouless and 
Landauer [!, 0] , it was possible to formulate the so- 
called one-parameter scaling theory of localization , one 
of the most fruitful approaches to the disorder-induced 
metal-insulator transition. The essential hypothesis of 
the scaling theory is that, close to the transition, a single 
relevant scaling variable describes the critical behavior. 
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An essential result of the one-parameter scaling theory 
is that the Anderson transition exists only in dimensions 
larger than two. In one dimension, all electronic states 
are localized, whatever the degree of randomness. In two 
dimensions, they are all localized, but only marginally, 
i.e. with a localization length exponentially large (thus 
possibly much larger than the sample size) for weak dis- 
order. 

In analogy to standard second-order phase-transitions, 
the localization length I is assumed to diverge at critical- 
ity according to a power law: 



i~ (w- w c y 



(2) 



with v the localization length critical exponent, W the 
disorder strength and W c the critical disorder strength. 
The most important assumption of the theory, the one- 
parameter scaling hypothesis, was numerically validated 
using a finite-size scaling method developed in 0, licj ]. 
This technique, which implements a real space renormal- 
ization, allowed to establish numerically the existence 
of a scaling function for the localization length. How- 
ever, the critical exponents measured using this method, 
v ~ 1.57 [ill . ll2T |. were not compatible with the result 
v = 1 obtained from a self-consistent approach of local- 
ization based on diagrammatic techniques, as developed 
in Ref. 

In the half-century since its birth, the Anderson model 
has become a paradigm for the studies of the interplay 
of quantum effects and disorder. Despite that, relatively 
few experimental results are available, for the following 
reasons: i) It is experimentally hard to finely tune the 
disorder in a real crystal; ii) the decoherence sources (col- 
lision with phonons, etc.) are difficult to master iii) 
electrons in a crystal present strong mutual interactions 
[3, (HI and iv) the wavefunction of the electrons in the 
crystal is not directly accessible, only transport proper- 
ties can be directly measured [la ]. 

It is thus interesting to search for other systems that 
display the Anderson transition, but are more favorable 
for experimental studies. Indeed, the concept of Ander- 
son localization has progressively been extended from its 
original solid-state physics scope to a variety of systems 
where a wave propagates in a disordered medium for ex- 
ample electromagnetic radiation [13, [HI and sound waves 
[TsL H3, HH . Photons propagating in disordered materials 
revealed to be an excellent system to observe the effects 
of localization However, in such systems, there is al- 
ways some absorption, whose signature can be quite sim- 
ilar to the signature of localization. Also, the measured 
quantity is the transmission, and the wavefunction itself 
is not accessible. The recent experimental observation of 
Anderson localization [12] using ultra-cold atomic matter 
waves has been done in a ID situation where states are 
always localized and no metal- insulator transition exists. 

A very interesting Anderson-type system is obtained 
by combining the Anderson model with another paradig- 
matic system, the kicked rotor (KR), which has been the- 
oretically studied for almost three decades. This system 



is well known to be classically chaotic [23|, and chaos 
plays here the role of a "dynamical" disorder. In the quan- 
tum case, the KR displays a localization phenomenon, 
called "dynamical localization" [24] which is analogous to 
the ID- Anderson localization [25J|. Moreover, a quasi- 
periodic generalization of the kicked rotor, substantially 
equivalent to the 3D Anderson model, was numerically 
shown to display an Anderson-like phase transition [2a ] . 
Experimental studies of the quantum kicked rotor were 
boosted by the realization of such a system with laser- 
cooled atoms interacting with a standing wave by Raizen 
and co-workers, which observed, for the first time, the 
Anderson localization with matter waves (27| . 

In the present paper we describe in detail a realiza- 
tion of an atomic matter-wave system that allows us 
to observe the Anderson metal-insulator transition [281 ]. 
We report a full characterization of this phase transition 
which includes an experimental validation of the one pa- 
rameter scaling hypothesis and the first non ambiguous 
experimental determination of the critical exponent v. 
Last but not least, we show numerically that the quan- 
tum chaotic system we consider has the same critical 
behavior as the true random 3D Anderson model. In 
particular, we show that the two models belong to the 
same universality class. Sec. [11 introduces the cold-atom 
realization of the periodic (standard) KR and its equiva- 
lence with the ID Anderson model, as well as the quasi- 
periodic generalization of this system that is equivalent 
of the 3D- Anderson model. Sec. IIIII describes the cor- 
responding experimental setup, paying attention to its 
experimental limits (decoherence, stray effects, limited 
observation time). In sec. HVI we report our direct exper- 
imental observation of the metal-insulator transition. In 
sec. |V] a scaling procedure is introduced that allows us 
to overcome experimental limitations and determine the 
critical exponent corresponding to the Anderson transi- 
tion. Sec. |Vl]is devoted to the universality of the critical 
behavior. Sec. IVIII concludes the paper. 



II. THE ATOMIC KICKED ROTOR AND ITS 
RELATION TO THE ANDERSON MODEL 

A. The atomic kicked rotor 

Consider a two level atom interacting with a laser 
standing wave of frequency ujl = k^c detuned by Al = 
ojl — ujq from the atomic transition of frequency ojq- It 
is well known that there are two kinds of interactions 
between the atom and the radiation: Firstly, the atom 
can absorb a photon from the laser and re-emit it spon- 
taneously in a random direction. This is a dissipative 
process giving rise to radiation pressure force, whose rate 
is ril 2 /4A| where T is the natural width and fl the res- 
onant Rabi frequency (we assume |Al| 3> T). Secondly, 
the atom can pick a photon in a laser mode and emit it in 
the same (or another) laser mode by stimulated emission. 
This conservative process is associated with a potential 



3 



acting on the atom's center of mass motion, called the 
optical or dipole potential. For a standing wave this po- 
tential is: 

Ml 2 

K P t = ^r- cos (2k L X) (3) 

where X is the atom center of mass position along the 
standing wave. Clearly, this interaction is one dimen- 
sional, as momentum exchanges between the atom and 
the radiation are always along the standing wave: The 
atom absorbs a photon in one of the propagating beams 
and emits it in the counterpropagating beam, leading 
to a quantized momentum exchange of 2fiJc£ along the 
X axis. An important point is that the optical poten- 
tial amplitude scales as 2 /Al whereas the spontaneous 
emission rate scales as rf2 2 /A|. In the regime |A/,| ^>r, 
the optical potential is the dominant contribution to the 
dynamics, with spontaneous emission events being rare. 
Moreover, one can reduce the spontaneous emission rate 
by increasing the detuning A^, provided that the laser 
has enough power to keep the potential amplitude at the 
required level. 

Suppose now that, instead of having the atom interact- 
ing continuously with the standing wave, one modulates 
the radiation intensity periodically (with period T\) so 
that it is on for a short time t (as compared to the atom 
dynamics) and off the rest of the period. One then ob- 
tains the Hamiltonian: 



k plays the crucial role of an effective Planck constant, 
which can be adjusted at will by modifying e.g. the pe- 
riod T\. As shown in the following, the most interest- 
ing physics takes place in the momentum. The scaling 
Eqs. |(5]) is such that P = 2hk]^ corresponds to p = k. 
If the atom is cold enough that its typical momentum is 
comparable to 2ftfc£ (the "quantum" of momentum ex- 
change), quantum effects can be observed in the system. 
Fortunately, magneto-optical traps produce atoms with 
a typical momentum of a few hk^. It is customary to 
measure the atomic momentum P in units of 2fifci, i.e. 
measure p in units of k. We thus will use: 

P P 

»=k = 2Hk- L ® 

For K > 5, the classical KR is fully chaotic, and the 
dynamics, although perfectly deterministic, behaves like 
a pseudo-random diffusive process known as "chaotic dif- 
fusion". For this reason, K is usually called "stochasticity 
parameter". The existence of classical chaos can be seen 
by integrating the classical equations of motion corre- 
sponding to Eq. © over a period, which leads to the 
so-called "Standard Map": 

xt+i —x t = pt (9) 
Pt+i-pt = A'sinxt+i. (10) 



p2 t,C)2 T ^ 

H = 2M + 8A7 COS {2kLX) £ ~ nT ^ ^ 

n 

where 5 T (t) = 1/t if \t\ < r/2 and zero otherwise. This 
functions tends to the Dirac (5-function as r — > 0. 

It is useful to introduce a set of scaled, dimensionless 
units: 



P 
t 

K 

k 

n 



2k L X 

2k L T 1 P/M 

t'/T, 

hn 2 T lT k 2 L 

2MA L 
M 



(5) 



In the limit of short pulses t<Ti, one then has: 



-w P 



K cos x S(t — n) 



(6) 



which is precisely the Hamiltonian of the kicked rotor 
[23L [2^l . One has thus realized an atomic kicked ro- 
tor [22]]. The above Hamiltonian is associated with the 
Schrodinger equation: 



ik 



dip 
~dt 



Hi). 



(7) 



If the stochasticity parameter K is large enough, sina^ 
generates random numbers for successive t values. The 
momentum then performs a random (though determin- 
istic) walk and the kinetic energy (averaged over the ini- 
tial conditions) increases linearly with time. If - as we 
assume in the following - the initial state is a narrow mo- 
mentum distribution centered around the origin p = 0, 
one obtains: 



(p 2 )(t) = Dt, 



(11) 



with D ?s K 2 /2k 2 being the diffusion constant. 

In the quantum chaotic diffusion is observed 

for times shorter than a characteristic "localization time" 
Tioc = D/2, after which quantum interferences build- 
up in the system that eventually "freeze" the dynam- 
ics, suppressing the diffusion. The mean kinetic en- 
ergy then tends to a constant (p 2 ){t — * 00) — > 2£ 2 with 
£ « K 2 /4k 2 At the same time, the momentum dis- 
tribution changes from a Gaussian shape characteristic 
of a diffusive process to a localized, exponential shape 
w exp(—\p\/£). This phenomenon is called "dynamical 
localization" (DL), "dynamical" meaning that the local- 
ization takes place in momentum space. In fact, as shown 
below, DL is intimately related to the Anderson localiza- 
tion, with, however, an important difference: DL takes 
place in momentum space, whereas Anderson localization 
is in real space. 
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B. Equivalence with the ID-Anderson model 

Let us consider the KR quantum dynamics. From a 
stroboscopic point of view, the motion is determined by 
the evolution operator over one period: 



U 



—iK cos x/h —ip /2k 



(12) 



whose eigenstates form a basis set allowing to calculate 
the temporal evolution. These Floquet states \<j>) are fully 
characterized by their quasienergy w, defined modulo 2ir: 



u\<t> u 



(13) 



The Hamiltonian, Eq. ([6]), is 27r-periodic in position x, 
and so is the evolution operator, Eq. (fl2|) . The Bloch 
theorem tells us that a Floquet eigenstate is a product of 
a periodic function of a; by a plane wave expiflx with < 
/3 < 1, is a constant, 0k being usually called the "quasi- 
momentum". A trivial transformation shows that one 
can equivalently consider periodic functions of x governed 
by the Hamiltonian, Eq. |(6|), where p is replaced by p+0k. 
In the following discussion, we will omit for simplicity the 
quasi-momentum, although it is straightforward to take 
it into account. Note that in all numerical simulations 
shown hereafter, we perform an averaging over the quasi- 
momentum, to follow the experimental conditions where 
an incoherent sum of all quasi-momenta is prepared. 

At this point, contact with a ID Anderson tight- 
binding model can be made by reformulating Eq. lfT3|) 
for the Floquet states j25|. Firstly, we rewrite the (uni- 
tary) kick operator: 



with 



-iK cobx/H = l+iW(x) 
l-iW{x) 



W{x) = tan(ifcosx/2S) . 



(14) 



(15) 



The periodic function W(x) can be expanded in a Fourier 
series: 



W(x)=J2 W re" 



Similarly, for the kinetic part, one gets: 



-i(p 2 /2k—uj) 



1 + iV 
1-iV ' 



(16) 



(17) 



with V diagonal in the momentum eigenbasis \m) = \p = 
km). Secondly, we make the following expansion in the 
momentum eigenbasis: 



(18) 



Then, the eigen-equation for the Floquet state can be 
rewritten: 

£m$m + W r^m-r = -Wb$fn , (19) 



with e m = tan - m 2 k/2)] [3l|. 

This is the equation for a tight-binding model with 
hopping elements W r to the r th neighbor, with eigen- 
energy Wo, and with on-site energy e m . The hopping el- 
ements are not restricted to nearest-neighbors, but they 
decrease exponentially with r [32]. In the original An- 
derson model, a random distribution is assigned to e m . 
Here, the sequence e m , although not satisfying the most 
stringent mathematical tests of randomness, is neverthe- 
less pseudo-random. These two conditions are sufficient 
for the Anderson localization to take place. The hopping 
integrals W r increase with the kick strength K, which 
thus plays the role of a control parameter in the Ander- 
son model fl9|) . Note that if k is a rational multiple of 
2tt, the e m are periodic in m. This leads to the quan- 
tum resonances of the kicked rotor, where the states are 
extended. 

When k is incommensurate with 2ir, the Floquet states 
are found to be exponentially localized, and this prop- 
erty accounts for dynamical localization. As shown 
in [33| , the localization length observed at long times for 
a wavepacket is essentially identical to the localization 
length of individual Floquet states. 

Many references discuss the detailed correspondence 
between quantum behavior of this dynamical system and 
Anderson localization: In Ref. [34[ an analogy between 
the KR and band random matrices was pointed out; the 
latter have been reduced to a ID nonlinear a model I35fl 
similar to those employed in the localization theory [3a |. 
In Ref. [13] the direct correspondence between the KR 
and the diffusive supersymmetric nonlinear a model was 
demonstrated. A diagrammatic approach [38] to the dy- 
namical localization in the Kicked Rotor was reported in 



C. The quasi-periodic Kicked Rotor and its 
analogy to the 3D-Anderson model 

As the Anderson transition exists only in three (or 
more) dimensions, one must generalize the KR to obtain 
a system analogous to a 3D Anderson model. 

Different generalizations of the KR have been the- 
oretically considered as analogs of the 3D-Anderson 
model [4oL [HI ]. Here we use the convenient three- 
incommensurate-frequencies generalization introduced in 
Refs. SHI: 



P 

H q p = — — h AT(t) cos x 2_] S(t 



(20) 



obtained simply by modulating the amplitude of the 
standing wave pulses with two new frequencies u>2 and 

K(t) = K [1 + e cos (u} 2 t + (p 2 ) cos (cj 3 t + (p 3 )} . (21) 

One can legitimately ask: where is the three dimensional 
aspect in the latter Hamiltonian? An answer can be given 
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by drawing a formal analogy between the quasiperiodic 
kicked rotor and a 3D kicked rotor with an initial condi- 
tion taken as "plane source" (see below). 

We start from the Hamiltonian of a 3D periodically 
kicked rotor: 

P 2 

H = — + UJ2P2 + U 3 p 3 

+K cos x\ [1 + e cos X2 cos £3] S(t — n) , (22) 

n 

let us consider the evolution of a wavefunction \& with 
the initial condition: 

$(xi,X2,X3,t = 0) = 5(0:1, t = 0)8(X2 - lf2)S(x 3 - 933) 

(23) 

The initial state being perfectly localized in X2 and 
x 3 , it is entirely delocalized in the conjugate momenta P2 



and p 3 , and can thus be seen as a "plane source" [3| in 
momentum space. 

From a stroboscopic point of view, the time-evolution 
of \& is determined by the evolution operator over one 
period: 

JJ _ g— iif cos a;i(l+£ cos x 2 008X3)/^ x g — i(pi/2+w 2 p2+^3P3)/K 

(24) 

It is then straightforward to see that the 3D-wave func- 
tion W at time t is related to its initial condition as: 

V(x 1 ,x 2 ,x 3 ,t) = U t V(x 1 ,x 2 ,x 3 ,t = 0) 
= 'B(x 1 ,t)5(x 2 -<pa- u> 2 t)5{x 3 - cp 3 - u> 3 t) , (25) 

with: 



S(xi t) = Y\ p - iK cosx i[ 1 + £cos ('P2+U2t')cos(ip 3 +tAj3t')]/k e -ipl/2ti i — Q) (26) 



On the other hand, consider now the evolution of an 
initial wave function ip(x,t = 0) with the Hamiltonian 
H qp of the quasiperiodic kicked rotor. It is also deter- 
mined by an evolution operator from kick to kick, but 
now this evolution operator U qp (t; t—1) depends on time, 
since the Hamiltonian H qp , Eq. (|20|) . is not time-periodic: 

U qp (t;t-1) = 

g — iK cos cos(cp2+^2^) cos(^3+W3i)]/Sg— ip 2 /2k (27) 

The wave-function ip(t) at time t is obtained by applying 
successively Uq P (t';t' — 1) for t' from 1 to t: 

t 

xjj(x, t)=]J U qp (t'; t' - l)^(x, t = 0). (28) 
t'=i 

From Eqs. ((25]), {27j) , J26]) and ((28]), it follows that 
tp{x,t) and H(xi,t) follow exactly the same evolution. 
Consequently, the dynamics of the quasiperiodic kicked 
rotor is strictly equivalent to that of a 3D kicked rotor 
with a plane source. Our experiment with the quasiperi- 
odic kicked rotor can be seen as a localization experi- 
ment in a 3D disordered system, where localization is 
actually observed in the direction perpendicular to the 
plane source [2l|. In other words, the situation is thus 
comparable to a transmission experiment where the sam- 
ple is illuminated by a plane wave and the exponential 
localization is only measured along the wave vector di- 
rection. Therefore, the behavior of the wave function ip 
subjected to the quasiperiodic kicked rotor Hamiltonian 
-ffqp, Eq. (|20j) . depicts all the properties of the dynamics 
of the quantum 3D kicked rotor, Eq. (|22|l . 



The Hamiltonian H, Eq. (|22|) . is invariant under the 
following transformation, product of time-reversal with 
parity: 

T : t -> -t,x -> -x,p -> p , (29) 

which is relevant for dynamical localization [H, 5i|. 
The evolution of the states according to the Hamilto- 
nian, Eq. (j20l) . is governed by the operator U, Eq. (]24p . 
which belongs to the Circular Orthogonal Ensemble class 
[13, SI, with the additional constraint at t = Eq. (|23| . 
Of course, the transformation (|29|l amounts to chang- 
ing (</J2,V3) to (— if2, — ^3) into the constraint {23]), i.e. 
to starting from a different wavefunction. On the other 
hand, from Eq. |25]l . one clearly sees that after t steps, 
the constraint reads S(x2 — ^2)8{x 3 — ip 3 ), with 

& = <P2 + U2t <p 3 = <fi3 + UJ 3 t. (30) 

Since the frequencies u>2 and ui 3 are incommensurate, the 
preceding equation immediately tells us that, along the 
time evolution, the constraint on the wavefunction can 
be arbitrary close to any phases (^2,^) [!§]■ This way, 
the time evolution results in an average over (almost) all 
possible phases, showing thus that the localization prop- 
erties are independent of a particular choice (cp2, ip 3 ), but 
only depend on the operator U. Therefore, the dynam- 
ical properties of the present quasiperiodic kicked rotor 
also belong to the orthogonal ensemble. 

It should be noted that the 3D aspect comes from the 
presence of 3 frequencies in the dynamical system: the 
usual "momentum frequency" k present in the standard 



kicked rotor Eq. ([6]) , and two additional time-frequencies 
UJ2 and cos. Thus, increasing the number of incommensu- 
rate frequencies allows one to tune the effective dimen- 
sionality of the system. 

Let us now consider the conditions for the observation 
of Anderson localization with the quasiperiodic kicked ro- 
tor. As for the standard kicked rotor, the Floquet states 
of the time-periodic 3D Hamiltonian H, Eq. |22"|) . can be 
mapped onto a 3D Anderson-like model: 



(31) 



where m = (mi, m%, 1713) and r label sites in a 3D cubic 
lattice, the on-site energy e m is 



{ - 

2 



. mi 



+ UJ 2 m 2 + W3TO3 



(32) 



and the hopping amplitudes W r are coefficients of a three- 
fold Fourier expansion of 

W(xi,X2, X3) = tan [K cos xi(l + e cos £2 cos £3) / 2Jc] . 

(33) 

An obvious necessary condition for the observation of 
localization effects is that e m is not periodic. This is 
achieved if (%, u)2,u>3, ir) are incommensurate. Of course, 
the presence of disorder in the diagonal energy e m is 
crucial to observe Anderson localization. When % is in- 
commensurate with 2ir, due to the presence of a non- 
linear dispersion in the mi direction, the classical dy- 
namics can become chaotic with diffusive spreading in 
all m directions [H, [H3|. A typical numerical simula- 
tion is shown in Fig. [TJ the classical motion is almost 
perfectly diffusive along the three Pi coordinates with a 
characteristic Gaussian shape in each direction. From 
Eq. |33l it is clear that hopping along the directions 
"2" and "3" is diminished by a factor e compared to 
hopping along direction "1". Not surprisingly, diffusion 
along pi is slightly faster than along p 2 and P3. The 
quasi-periodically kicked rotor is thus analogous to an 
anisotropic Anderson model [5ll . I52L [53| . 

When those conditions are verified, localization effects 
as predicted for the 3D Anderson model are expected, 
namely either a diffusive or a localized regime. Localized 
states would be observed if the disorder strength is large 
as compared to the hopping. In the case of the model 
Eq. (|3"Tjl . the amplitude of the disorder is fixed, but the 
hopping amplitudes can be controlled by changing the 
stochasticity parameter K (and/or the modulation am- 
plitude e): W r is easily seen to increase with K. In other 
words, the larger K, the smaller the disorder. One thus 
expects to observe diffusive regime for large stochastic- 
ity or/and modulation amplitude (small disorder) and 
localized regime for small K or/and e (large disorder). 
It should be emphasized that there is no stricto sensu 
mobility edge in our system. Depending on the values of 
the parameters K,k,£,u>2,k>3, all Floquet states are lo- 
calized or all are delocalized. The boundary of the metal- 
insulator transition is in the K, H, e, o; 2 , w 3 space. As seen 
below, K and e are the primarily important parameters. 
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Figure 1: (Color on line) Classical diffusive motion for the 3D 
kicked rotor Eq. 1221 The initial state is localized around the 
origin. After 1000 kicks, the classical momentum distribution 
(black and red curves) has the Gaussian shape characteris- 
tic of a diffusive motion. The blue and green curves are fits 
by a Gaussian which do not show any statistically significant 
deviation. The black (resp. red) curve is the momentum 
distribution along pi (resp. pa). The distribution along p3 
is identical that along P2- The anisotropic diffusion happens 
because the hopping along the directions "2" and "3" is di- 
minished by a factor e compared to hopping along direction 
"1". Parameters are K = 10, J5 = 2.85, e = 0.8, LJ 2 /27r = y/E, 
lu 3 /2it = a/13. 



In an analytical work on a similar problem [43] Basko et 
al. showed that the weak dynamical localization regime 
of a d-frequency quantum dot system is similar to the 
weak localization in a d-dimensional Anderson model. 
This work confirms the equivalence between our system 
and the 3D- Anderson model. The above arguments were 
also validated numerically (42L ]5p| . 

Numerical simulations of the evolution of the quasi- 
periodically kicked rotor are straightforward. The free 
evolution between consecutive pulses is diagonal in mo- 
mentum representation, while the kick operator is di- 
agonal in position representation (whatever the kick 
strength, constant or quasi-periodic). Switching between 
momentum and position representation is easily done 
through a Fast Fourier Transform. We are thus able to 
compute the evolution of a large number of initial states 
(typically one thousand) over a very long time (typically 
up to one million kicks, much more than in the exper- 
iment). All numerical results shown below have been 
carefully checked for convergence. Except when explic- 
itly stated, averaging over the quasi-momentum (3 has 
been performed, in accordance with the experimental re- 
alization. 



III. EXPERIMENTAL REALIZATION WITH 
ATOMIC MATTER- WAVES 

A. Experimental setup 

Our experimental setup has been described in detail 
in previous publications [Hil . [5ol [56l . [H3] and was used in 
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various inv estig ations on the quasiperiodic kicked rotor 
[Hsl . [59l . 16(1 l6ll f Briefly, our experiments are performed 
with cesium atoms produced in a standard magneto- 
optical trap (MOT). A long Sisyphus- molasses phase (25 
ms) allows us to obtain 10 7 atoms at a measured temper- 
ature of 3.2 fiK. The velocity distribution of the atoms is 
well modeled by an incoherent sum of plane waves form- 
ing a Gaussian of full width at half maximum (FWHM) 
equal to 8hkL, which is much narrower than the ex- 
pected localization length. The MOT beams and mag- 
netic field are turned off and the sequence of kicks is 
applied to the atoms. The beam forming the standing 
wave passes through an acousto-optical modulator driven 
by RF pulse synthesizers, which generates the kicks at a 
typical frequency of 1/Ti = 36 kHz (which corresponds 
to % = 2.89), of duration r =900 ns and with a raising 
time of 50 ns. The modulation is thus an almost perfect 
square, at the time scale of the atomic motion, and its 
duration and period can be set by a microcomputer. The 
beam is then injected in an optical fiber that brings it to 
the interaction region, and the standing wave is obtained 
simply by back-reflection of this beam. The standing 
wave has a typical power 160 mW, its profile intensity 
has a FWHM of 1.5 mm, and it is far off-resonant (7.3 
GHz to red, or 1.4 x 10 3 r), in order to reduce spontaneous 
emission. The corresponding stochasticity parameter is 
K ps 15. 

A very interesting property of our system (as com- 
pared to solid-state systems) is that the wave function 
is accessible (or at least its square modulus). We mea- 
sure the atomic velocity distribution by velocity-selective 
Raman stimulated transitions, which are sensitive to the 
atomic velocity via Doppler effect, allowing an optimal 
velocity resolution of about 2 mm/s. A Raman pulse de- 
tuned of Sr with respect to the Raman resonance trans- 
fers the atoms in the velocity class v — <5j?/ (2fcjj) — vr 
with vr = hkn/M (Jcr is the wave number of the Ra- 
man beams) from the F g = 4 to the F g = 3 ground- 
state hyperfine sublevel. A beam resonant with the 
F g = 4 — > F e = 5 transition is then applied to push 
the remaining atoms out of the interaction region. The 
F g = 3 atoms are then optically pumped to the F g = 4 
sublevel and interact with a resonant probe beam: The 
absorption signal is thus proportional to the population 
of the F g = 4 level, thus to the population of the selected 
velocity class. The whole sequence then starts again with 
a different value of the Raman detuning to probe a new 
velocity class, allowing a reconstruction of the velocity 
distribution [5H.IB7I]. 



B. Decoherence sources 

Any quantum experiment must consider decoherence 
sources that destroy quantum interference effects (in our 
case, localization) reestablishing a diffusive dynamics. 
The most important sources of decoherence in our exper- 
iment are (i) atomic collisions, (ii) spontaneous emission, 



and (iii) the deviation of the standing wave from strict 
horizontality. 

For an isolated system described by a single wavefunc- 
tion, phase coherence between different positions is "per- 
fect". When the system is weakly coupled to an exter- 
nal bath, it cannot be any longer described by a single 
wavefunction; the most convenient description usually in- 
volves a density matrix p. Non-diagonal matrix elements 
of the type (x\p\x') quantify the degree of coherence of 
the system between position x and x' . As a general rule, 
the effect of the external bath is to make the non-diagonal 
elements of the density matrix to decay relatively rapidly, 
more rapidly than the diagonal elements: this is decoher- 
ence (not to be confused with dissipation) jfj^]. Effects 
like Anderson localization are due to subtle destructive 
interference amongst various components of the wave- 
function, which inhibit the classically allowed transport: 
they are thus very sensitive to decoherence. One usually 
quantify the strength of decoherence effects by defining a 
phase coherence time, the characteristic time over which 
the non-diagonal elements of the density matrix decay be- 
cause of coupling to the external bath. In our case, the 
non-diagonal element of interest are between eigenstates 
\p) and \p') located at a typical distance \p — p'\ compa- 
rable to the localization length in momentum space. 

Localization effects can be observed only for times 
shorter than the phase coherence time Beyond the 
phase coherence time, interference effects are killed and 
classical-like diffusive dynamics sets in. In the following, 
we shall express the characteristic times of the decoher- 
ence processes (i), (ii) and (iii), as functions of the ex- 
perimental parameters to show that they can be set large 
enough for localization effects to be observable. 

In atom-atom collisions, the dominant effect is that of 
collisions between cold atoms, the density of the cloud 
being around 8 orders of magnitude larger than the den- 
sity of the background hot gas. A cloud density of 10 12 
cm~ 3 with a mean velocity 1 cm/s and a collision cross- 
section of 6 x 10 -11 cm 2 gives a collision rate of ps 60 
s" 1 , or 1.6 x 10~ 3 per kick; the collision phase coherence 
time is thus ~ 600 kicks. 

In order to have a better idea of the decoherence ef- 
fect induced by spontaneous emission, let us consider 
the temporal evolution of an initial plane-wave function: 
ip(p,t = 0) = S(p — p ) evolving with the KR Hamilto- 
nian Eq. ([6]). After dynamical localization sets in, the 
momentum distribution ceases to expand because of de- 
structive interference between the various components of 
the wavefunction. Spontaneous emission brings a random 
recoil to the atomic momentum which is not an integer 
multiple oi2Kki J . Thus, the quasi- momentum j3 performs 
a random jump. As the phase factors involved in the free 
evolution depend on the quasi-momentum, the relative 
phases between interfering paths are scrambled, result- 
ing in a new transient diffusive behavior for another du- 
ration of t 1oc . DL is thus expected to be destroyed if 
spontaneous emission is regularly repeated. Note that a 
single spontaneous emission event completely breaks the 
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Figure 2: (Color online) Gravity effects on a slightly inclined 
kicked rotor Eq. l134[| . The deviation of the standing wave 
from horizontality is a = 0° (black lower curve), a = 0.1° 
(red lower middle curve), a = 0.4° (green upper middle curve) 
and a — 1° (blue upper curve). The stochasticity parameter is 
taken as K = 5 and the effective Planck constant is £ = 2.85. 
The dynamics of an initial thermal state is simulated and the 
corresponding mean kinetic energy is plotted versus time. For 
angles larger than 0.1°, the slow drift of momentum induces 
a diffusive behavior clearly visible on the time-scale of the 
experiment. 



drift of the quasi- momentum at constant rate —rig, whose 
effect is to break dynamical localization. Indeed, the 
destructive interference between various components of 
the momentum wavefunction - responsible for dynam- 
ical localization - is partially destroyed by the quasi- 
momentum drift, as the various phase factors accumu- 
lated during the free evolution between two consecutive 
kicks, exp [— i(m + /3S) 2 /2S] also drift. The net result is 
a residual diffusion constant, depending on rj g . Although 
this is not strictly a decoherence effect (the whole evo- 
lution is fully phase coherent), it similarly destroys dy- 
namical localization. We thus define the phase coherence 
time T g as the time needed to double (p 2 ) compared to 
the dynamically localized situation. Numerical simula- 
tions taking into account the gravity effect confirm the 
discussion above, see Fig. [2J If the standing wave devi- 
ates from horizontality by an angle a = 1°, then r g m 120 
kicks whereas when the angle a = 0.1°, r g « 350 kicks. 
In the timescale of the experiment (150 kicks), the de- 
viation from horizontality must be less than 0.1°. This 
decoherence effect is rather important. To the best of 
our knowledge, its importance was not fully appreciated 
in previous experiments. A detailed discussion of this 
effect will be presented elsewhere [64| . 



phase coherence, implying that the phase coherence time 
is simply the inverse of the spontaneous emission rate. 

Spontaneous emission tends to reestablish a diffusive 
evolution with a diffusion constant that is roughly Vjk 2 
where rj — T£1 2 t /&A 2 L , is the spontaneous emission rate 
expressed in photons per kick, which can be cast in the 
more useful form rj = (IY/8) (I/I s ) (T/Al) , where / 
is the intensity and I s « 2.2 mW/cm 2 is the transition 
saturation intensity. Around the transition (K ~ 6), the 
experimental values indicated above give rj « 2.1 x 10~ 3 
s , or a typical phase coherence time of ~500 kicks. 

Another effect leading to the destruction of localization 
is the standing wave deviation from horizontality. In this 
case, a gravity term must be added in the Hamiltonian 

H g = y - r) g x + K cos xy]6(t- n) , (34) 

n 

The dimensionless gravity term r) g is: 

V ° = 2hki n8ma ' (35) 

with g the gravity acceleration and a the angle between 
the horizontal direction and the standing wave. The 
physical interpretation is quite clear: mgTisina is the 
additional momentum transferred to the atoms between 
two consecutive kicks, which must be compared to the 
width of the Brillouin zone 2ftfci. 

The gravity term — r) g x breaks the spatial periodic- 
ity of the Hamiltonian, and consequently the conserva- 
tion of the quasi- momentum (3k. It actually produces a 



C. Conditions for the observation of localization 
effects 

We now discuss the conditions that must be satisfied 
in order to observe localization effects experimentally. 

Firstly, the system must present some kind of disorder: 
As discussed in section III C| this means that te, u>2 and 
W3 and 7r must be incommensurate. This is achieved if 
we take k = 2.89, uj 2 = and W3 = 2iry/l3. A 

more detailed discussion concerning the choice of these 
parameters will be given in section [VI Bl 

Secondly, in order to observe dynamical localization 
effects instead of trivial classical localization, we must be 
in a regime where the classical system has no KAM bar- 
riers which can prevent the classical diffusive transport. 
For the standard, periodic KR, full chaos is obtained for 
K > 4. In order to determine the corresponding thresh- 
old for the quasiperiodic system, we performed numer- 
ical simulations of the classical dynamics corresponding 
to Eq. (J22j) , for various values of the stochasticity pa- 
rameter. The dynamics is found to be fully diffusive for 
K > 2, a considerably smaller value than for the stan- 
dard KR. In particular, no classical localization effects 
due to KAM barriers are observed for K > 2. In any 
case, the experiments and the numerical simulations in 
the following are all performed for K > 4, where the 
classical dynamics is diffusive, see Fig. [31 

Thirdly, short enough pulses must be used that they 
can be considered as delta pulses (65|. Numerical sim- 
ulations of the quasiperiodic kicked rotor with a finite 
pulse duration r = 0.9/is and a thermal initial momen- 
tum distribution show that less than 1% of the atoms are 




Figure 3: Classical chaotic diffusion for the quasiperiodic 
kicked rotor Eq. lf22]) . The dynamics of an initial thermal 
distribution of classical particles is simulated and the corre- 
sponding mean kinetic energy is plotted versus time (number 
of kicks). The stochasticity parameter K (the modulation 
amplitude e) varies linearly between 4 and 9 (0.1 and 0.8), 
following the experimental path, Fig. [6] The dashed line of 
slope 1 demonstrates the linear increase of (p 2 ) vs. time t. No 
classical localization effects are observed. The chaotic diffu- 
sion is characteristic of the presence of pseudo-disorder in the 
quasiperiodic kicked rotor, leading to a pseudo-random walk 
in momentum space. 




sensitive to the duration of the pulses. Only atoms in 
the tails of the momentum distribution have sufficiently 
large atomic velocity to move by a significant fraction of 
Xl during the pulse, thus feeling a smaller effective kick. 

Fourthly, a sufficiently narrow initial state must be pre- 
pared in order to observe dynamical localization, i.e. the 
freezing of the initial diffusive expansion of the wave- 
function into an exponentially localized state. A suffi- 
cient condition is that the initial width of the momentum 
distribution be smaller than the localization length. In 
our system, we have an initial momentum distribution of 
half- width 2k. This is comparable to the shortest localiza- 
tion length at the lowest K = 4 value, as experimentally 
proved, see inset of Fig.[5l A consequence is that, in this 
regime, the exponential shape of the wavefunction after 
dynamical localization is established is slightly rounded 
at the tip. For higher values - say K > 5, - the initial 
width of the atomic wavefunction can be safely neglected. 



Finally, decoherence processes must be kept small dur- 
ing the experiment. The large detuning of the standing 
wave allows to keep the spontaneous emission rate very 
small, i.e. the corresponding phase-coherence large as 
compared to the duration of the experiment. A good 
control on the horizontality of the standing wave insures 
that gravity do not lead to a destruction of localization 
effects on the time-scale of the experiment. 



Figure 4: (Color online) Experimentally measured momentum 
distributions after 150 kicks, exponentially localized in the 
insulator region (blue) and Gaussian in the diffusive (metallic) 
region (red), (a) linear scale, (b) log scale. For both curves 
k = 2.89, for the localized distribution (blue) K = 5.0 and 
e = 0.24, for the Gaussian distribution (red) K = 9.0 and 
e = 0.8. 



IV. EXPERIMENTAL OBSERVATION OF THE 
DISORDER INDUCED METAL-INSULATOR 
TRANSITION 



In a typical experimental run, we apply a sequence of 
kicks to the atomic cloud and measure its dynamics. In 
the localized regime, the evolution of its momentum dis- 
tribution is "frozen" after the localization time (typically 
of the order of 12 kicks at low K) into an exponential 
curve exp(— |p|/£). In the diffusive regime, the initial 
Gaussian shape is preserved and the distribution gets 
broader as kicks are applied, corresponding to a linear 
increase of the average kinetic energy. Figured] shows the 
experimentally observed momentum distributions, an ex- 
ponentially localized distribution for small K and e (blue 
curve), characteristic of dynamical localization, and a 
broad, Gaussian-shaped distribution for large K and e 
(red curve), characteristic of the diffusive regime. 

Measuring the whole momentum distribution takes too 
much time: one must repeat the whole sequence (from 
the preparation of a new atom cloud up to the Raman 
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Figure 5: (Color online) Temporal dynamics of the quasi- 
periodic kicked rotor. We experimentally measure the popu- 
lation Ilo(£) of the zero-momentum class as a function of time 
(number of kicks) and plot the quantity II^" 2 (t) oc (p 2 )(i). 
Clearly, it tends to a constant in the localized regime (blue 
lower curve corresponding to K = 4 and e = 0.1) and in- 
creases linearly with time in the diffusive regime (red upper 
curve corresponding to K = 9 and e = 0.8). The inset shows 
the behavior close to the localization time, k = 2.89. 

measurement of the velocity distribution) for each ve- 
locity class. Moreover, for each time step, a complete 
momentum distribution must be measured. Fortunately, 
it is sufficient, and much easier, to measure the popula- 
tion Ilo(i) of the zero velocity class, as IIq 2 {t) is propor- 
tional to (p 2 ){t) (the total number of atoms is constant). 
The proportionality factor between ll^ 2 (t) and (p 2 )(£) 
depends on the detailed shape of the momentum distri- 
bution and is thus different in the localized and diffu- 
sive regime, but this small difference is a small correction 
to the main phenomenon: divergence of the localization 
length near the transition. 

Note that, strictly speaking, the proportionality be- 
tween HQ 2 (t) and (p 2 )(t) breaks at criticality due to the 
multifractal character of critical states jsH]. However, on 
the time scale of the experiment (£ = 150 kicks), the de- 
viation from strict proportionality is seen (numerically) 
to be negligible. At longer times (thousands or millions 
of kicks) , the effect of multifractality is visible and quan- 
titatively measurable. This is beyond the scope of this 
paper and will be analyzed elsewhere (67| . 

For each run, a value of Ho(t) is recorded after a given 
number of kicks is applied, then the measurement se- 
quence starts again with the next number of kicks. We 
also record the background signal obtained by not apply- 
ing the Raman detection sequence, and the total number 
of atoms in the cold-atom cloud. These signals are used 
to correct the experimental data from background signals 
and long-term drifts of the cloud population. 

Figure [5] shows the experimentally measured IIq 2 {t) in 
the localized and diffusive regimes. It clearly shows the 
initial diffusive phase and the freezing of the quantum 
dynamics in the localized regime (blue curve in Fig. [5]) . 
Along with the observation of an exponential localization 
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Figure 6: (Color online) Phase diagram of the quasiperiodic 
kicked rotor, from numerical simulations. The localized (insu- 
lator) region is shown in blue, the diffusive (metallic) region is 
shown in red. The experimental parameters are swept along 
the diagonal dash-dotted line. 

of the wave-function in Fig.JU this constitutes a clear-cut 
proof of the observation of dynamical localization. In the 
diffusive regime, IIq 2 {t) is seen to increase linearly with 
time (red curve in Fig. [5]), corresponding to the Gaussian 
red curve in Fig. [H 

After having observed Anderson localization for strong 
effective disorder strength and diffusive transport for 
small effective disorder, the next step is to walk the way 
between these two regimes, and explore the phase transi- 
tion expected (numerically) to take place along a critical 
line in the plane (K, e > 0) (Fig. [6]) . In order to con- 
fine the transition to a narrow range of parameters, we 
choose a path that cross the critical curve (Fig. [6]) "at a 
right angle"; we thus vary simultaneously K and e along 
a line going from K = 4, e = 0.1 in the localized region 
to K = 9, e = 0.8 in the diffusive region; the critical line 
is then crossed at K = K c = 6.6. 

A simple way to investigate the phase transition is 
the following [42]. In the localized regime, wait for a 
time longer than the localization time so that a local- 
ized frozen wave-function is observed, then measure its 
localization length. One can in such a way study the 
behavior of the localization length vs. disorder: at criti- 
cality, it should diverge as I ~ (K — K c )~ v . This would 
give the critical stochasticity parameter K c and the criti- 
cal exponent v. However, we cannot proceed that way in 
our case, because when one approaches the critical point 
from the insulator side, the localization time diverges as 
n ~ £ 3 r~j (K — K c )~ 31 ' in three dimensions (see below). 
In our system, a localized momentum distribution would 
be observable in the vicinity of the transition only for 
prohibitively large numbers of kicks, which are, in prac- 
tice, limited to 150, essentially because of decoherence 
effects and because the free fall of the atom cloud takes 
it out of the standing wave. Consequently, it is vain to 
investigate experimentally the Anderson transition only 
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from static properties such as the divergence of the lo- 
calization length at criticality, which could be obtained 
only for t S> n oc . Fortunately, there is another way to 
observe the Anderson transition, which we shall present 
in the following sections. 



V. CHARACTERIZATION OF THE 
ANDERSON PHASE-TRANSITION 

Finite-time effects act as finite-size effects do on finite- 
size samples subjected to phase-transitions. Numerical 
simulations of the Anderson transition on the standard 
3D- Anderson model are necessarily performed on finite- 
size samples of finite size L. In the vicinity of the transi- 
tion, the localization length £ [see Eq. ([2])] diverges and 
thus can greatly exceed L. In this regime, L acts as an up- 
per bound for the effectively observed localization length 
th- This smooths the transition, no divergence of the lo- 
calization length can be directly observed on a finite-size 
sample. In order to overcome this limitation, a power- 
ful real-space renormalization method, called finite-size 
scaling [3, LLIJ], was introduced. This method is based on 
a single parameter scaling hypothesis 0| and allows to 
extrapolate from the scaling behavior oi II versus L the 
asymptotic value of the localization length I correspond- 
ing to L — > co. We can generalize static scaling laws to 
cover our time-dependent problem (see [68j] for a similar 
approach in percolation theory). The single parameter 
scaling theory [§], successfully used for the standard 3D 
Anderson model @, [lfj, can be applied to analyze our 
experimental and numerical data, and especially to de- 
termine the critical properties of the Anderson transition 
that we observe, i.e. the critical exponents. 



A. Scaling law at finite time 

Knowing the asymptotic behavior when t — ► oo is 
not enough, an additional time-dependent property is 
needed, too, which we shall investigate now. For K 
far above K c one observes normal diffusion, (p 2 ) oc t, 
whereas for K far below K c , the quantum dynamics 
freezes, at sufficiently long times. Following the stan- 
dard analysis of the Anderson transition, we make the 
hypothesis that the transition that we observe for the 
quasi-periodically kicked rotor follows a one-parameter 
scaling law [6S| (the validity of this scaling hypothesis will 
of course be checked at the end of the analysis) . At the 
critical point, a third kind of dynamics, namely anoma- 
lous diffusion, with (p 2 ) ~ t k k ^ 1, is expected. Let us 
consider in greater detail the behavior very close to K c 
where these three different laws merge. 

In the localized regime, for sufficiently long times, the 
behavior depends only on the localization length which 
diverges as K goes to K c : 

(p 2 ) ~ e ~ {K c - K)~ 2v (for K<K C ), (36) 



with v the localization length critical exponent. 

For K > K c , the mean kinetic energy increases lin- 
early with time, and the proportionality constant is the 
diffusion coefficient D(K). For K < K c , (p 2 ) is bounded 
by Eq. ((361 and there is no diffusion. Thus D(K) van- 
ishes below K c . A different critical exponent s is used to 
describe how D(K) goes to zero above threshold: 

D(K) ~{K- K c ) s (for K > K c ). (37) 

We shall now find a single expression presenting these 
two limit behaviors and also displaying anomalous diffu- 
sion at the critical point. We note that, according to the 
theory of phase-transitions in finite-size samples, a scal- 
ing can be applied to (p 2 ) depending on the two variables 
1/t and (K — K c ), both going to zero. We thus use the 
general scaling law: 

(p 2 ) = t^F[{K-K c )t k *] , (38) 

with F(x) an unknown scaling function. The exponents 
ki and k 2 can be determined as follows. 

In the diffusive regime, for long enough times, we must 
recover the diffusion law with D ~ (K - K c ) s [Eq. (|37|) ]: 
hence, for x 3> 1, the scaling function F(x) should scale 
as x s : 

As in the diffusive regime, (p 2 ) ~ t, we must have k\ + 
sk 2 = 1. 

In the localized regime, on the other hand, one must 
recover (p 2 ) - (K c - Ky 2 " [Eq. (J36j)] for sufficiently long 
times. Thus, for x — > — oo, F(x) — > (— x)~ 2v , giving: 

(p 2 ) = t fc i- 2 ^ (K c - K)- 2v (40) 

which is compatible with Eq. lj36|l only if ki = 2vk 2 - 
These two relations determine k\ and in terms of the 
physically more meaningful critical exponents s and v. 



In the standard Anderson model, the critical exponents 
are related by Wegner's scaling law @: 

s = {d-2)v, (41) 

with d being the dimensionality of the system. For our 
system, one obtains: 

fci = 2/3; k 2 = l/3u. (42) 

We therefore expect at the critical point anomalous diffu- 
sion with (p 2 ) = t kl F{0) - i 2 / 3 . We present in the next 
sub-section a numerical and experimental validation of 
this prediction. 
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Figure 7: (Color online) Numerically simulated time- 
evolution of (p 2 ) for the quasiperiodic kicked rotor. At the 
critical point K = K c « 6.4 (purple middle curve), anomalous 
diffusion (p 2 ) ~ i 2//3 is clearly observed, as expected from the- 
oretical arguments (cf. text). The log-log plot of the critical 
curve is very well fitted by a straight line of slope 0.664 (black 
dashed line). In the vicinity of the transition, the dynamics 
departs from the anomalous diffusion to tend gradually ei- 
ther to a diffusive dynamics (red upper curves corresponding 
to K => K c bending upwards for large t) or to a localized 
dynamics (blue lower curves corresponding to K < K c bend- 
ing downwards for large t). Other parameters are k = 2.85, 
u>2 = 2-7r\/5 and ljs = 2-7r\/l3. 
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B. Critical anomalous diffusion 

We verified numerically that the critical behavior, cor- 
responding to the anomalous diffusion in t 2 ' 3 is observed 
up to a very large number of kicks (t = 10 6 ). The (pur- 
ple) middle curve of Fig. [7] displays the time-evolution 
of (p 2 ) from numerical simulations for the stochasticity 
parameter K = 6.4. Anomalous diffusion (p 2 ) ~ t 2 l z is 
clearly seen from the log-log plot over 4 orders of magni- 
tude, which is very well fitted by a straight line of slope 
0.664. Other curves, for different K, tend at long times 
to bend either horizontally (below K c ) or towards slope 
unity (above K c ). This is a clear proof that we face here a 
true phase transition and not a smooth cross-over. Note 
also that the fact that the numerically measured critical 
slope is very close to the theoretical prediction 2/3 im- 
plies that the Wegner's scaling law s — v is valid at an 
accuracy better than 1%. 

Fig. [8] displays the experimental evolution of n^~ 2 (t) ~ 
(p 2 ) versus time. The critical curve (middle curve corre- 
sponding to K sa 6.4) in purple is well fitted by the rela- 
tion HV 2 (i) = A + Bt 2/3 , see Fig. Hi. Fig. [8b displays in 
log- log scale the experimental data Ilg 2 {t) vs t. The alge- 
braic dependence (with exponent 2/3) of the critical dy- 
namics is again clearly visible. In all plots in Figs. [8] the 
red upper curves evidence the above-criticality diffusive 
behavior and the blue lower curves the below-criticality 
localized behavior. 

From renormalization theory, we know that the crit- 



Figure 8: (Color online) Experimentally observed time- 
evolution of II^ 2 ~ <^p 2 ) for the quasiperiodic kicked ro- 
tor. Close to the critical point K = K c w 6.4 (purple mid- 
dle curve), anomalous diffusion il ( 7 2 (t) ~ t 2 ^ 3 is clearly ob- 
served, (a) The critical anomalous curve is well fitted by 
IIo 2 {t) =A + Bt 2/3 (black dashed line). The red upper curve 
evidence the far-above-criticality diffusive behavior (K = 9.0) 
and the blue lower curve the far-below-criticality (K = 4.0) 
localized behavior, (b) These experimental results show a 
clear algebraic behavior, with exponent ~ (blue lower curve, 
localized regime), 2/3 (purple middle curve, critical regime) 
and 1 (red upper curve, diffusive regime), slightly perturbed 
by decoherence processes responsible for the residual increase 
in the localized regime. Other parameters are the same as in 
Fig. El 



ical behavior shows the existence of a fixed hyperbolic 
point [16]. It is a fixed point because the critical behav- 
ior remains the same at all times (opposite to the local- 
ized case for example, for which a characteristic time can 
be defined, the localization time), and it is a hyperbolic 
point since the localized dynamics close to criticality will 
follow only for a finite time the anomalous diffusion with 
exponent 2/3 and will progressively tend to a localized 
behavior for large enough time. The rate at which the 
behavior changes is related to the critical exponent of the 
phase transition v. 

An efficient way to observe the departing of the dynam- 
ics from the critical anomalous diffusion is to consider the 
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Figure 9: (Color online) Numerical simulation showing the 
evolution of the dynamics from the critical behavior towards 
either a diffusive dynamics or a localized state. Plotting the 
quantity In A = ln((p 2 }t~ 2//3 ) vs. hit allows to easily distin- 
guish the critical behavior from diffusive or localized behavior: 
The critical curve (corresponding to K = Kc « 6.4) has a 
zero slope; whereas the far localized (K — 4.0) one has a slope 
—2/3 and the far diffusive (K = 9.0) one a slope 1/3. Other 
parameters are the same as in Fig. \7\ 




-5 5 10 

ln(1/t 1/3 ) lnfc/t 1/3 ) 

Figure 10: (Color online) Raw numerical data, displayed in 
the form In A = ln((p 2 )t _2/3 ) vs. lnt _1/3 (on the left). Each 
curve corresponds to a different stochasticity parameter K. 
The finite-time scaling procedure consists in shifting horizon- 
tally each curve by a quantity ln£(A") so that the curves over- 
lap. This allows one to determine both the scaling function / 
(on the right) and the scaling parameter £(K). 



quantity 



A 



<p!> 

t 2/3 



(43) 



or, equivalently, in the case of the experimental data: 



A, : 



1 



ng(t)tV3 



(44) 



as a function of time. This is illustrated in Fig. 0J which 
displays In A vs. t. The critical behavior can be eas- 
ily pin-pointed: The corresponding (purple) curve has a 
zero slope, as the quantity A is (asymptotically) constant 
at criticality. In the diffusive regime, the quantity A in- 
creases with time (red curves), whereas it decreases in 
the localized regime (blue curves). In the far localized 
regime, we observe an algebraic dependence A(i) ~ i~ 2 / 3 
as (p 2 ){t) = 2£ 2 for t > t\ oc . In the far diffusive regime, 
the algebraic dependence is A(t) ~ t 1 / 3 as (p 2 )(t) ~ t. 

The above numerical and experimental observations 
validate the theoretical prediction for the critical behav- 
ior: (p 2 ) ~ t 2 / z . Such critical behavior for the quasi- 
periodic kicked rotor was predicted using another scaling 
approach and numerically verified in [50j]. It was also 
numerically observed for a spatially-3D kicked rotor [!l| 
and in the standard 3D Anderson model jfjjj], and put on 
firm theoretical grounds in (70| . 



C. Finite-time scaling 

We shall now explain the procedure used to verify the 
scaling of our numerical and experimental data according 
to the law deduced above: 



- +2/3 



(K - K c )t 



(45) 



Our method is similar to the finite-size scaling procedure 
used by MacKinnon and Kramer 0, [7l[ , and Pichard and 
Sarma to numerically study the Anderson transition 
on finite-samples of the 3D Anderson model, but we ap- 
ply it here to the temporal behavior of the data, thus the 
name "finite- time scaling". 

We assume the quantity A(K,t) = (p 2 )i~ 2 / 3 to be an 
arbitrary function 



A(K,t) = f (aK)tr 1/3 ) , 



(46) 



where the scaling parameter £,(K) depends only on K, 
which is the parameter appearing in the one-parameter 
scaling hypothesis. This scaling assumption is less re- 
strictive than Eq. f45|) since no assumption on the de- 
pendence of £ on K is made. We must thus show that 
the resulting scaling parameter £(K) is compatible with 
Eq. (45]). 

In the left part of Fig. [10] we display plots of In A(K, t) 
vs. lnt^ 1 / 3 for different values of K. For most values of 
In A, several values of lnt -1 / 3 correspond to the same K 
value. The only way to conform with the condition f46]) 
is to shift each curve horizontally by a different quantity 
In £ (K) such that curves corresponding to different values 
of K overlap. This can be achieved by minimizing the 
variance of the values ln£(AT)t -1 / 3 corresponding to each 
value of In A. The function £(K) can be determined by 
applying a least square fit to the data. 

This minimization procedure does not allow one to 
compute the absolute scale of £{K), as the shifting proce- 
dure (see Fig. [TO]) is invariant under a global shift of the 
origin. We can thus set the scaling parameter £(K) to 
be equal to the localization length in the strongly local- 
ized regime where the duration of the experiment is much 
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Figure 11: (Color online) Finite-time scaling applied to the 
results of numerical simulations of the quasiperiodic kicked 
rotor. The time-evolution of (p 2 ) is computed as a function 
of time, from 30 to 10 4 kicks, for several values of K between 
K = 4 and K = 9. The finite-time scaling procedure al- 
lows us to determine both the scaling function / (a), clearly 
displaying an upper branch (red) associated with the diffu- 
sive regime, and a lower branch (blue) associated with the 
localized regime. The dependence of the scaling parameter 
£ on K (b) displays a divergent behavior around the critical 
point K c = 6.4, which is the signature of the Anderson phase 
transition. The dashed line is a fit using Eq. Q48|l. The re- 
sulting critical exponent is v = 1.6 ± 0.1. Other parameters 
are k = 2.85, u>2 = 1-k\J1> and W3 = 27r^l3. 



larger than the localization time, and (p 2 ) converges to 
its asymptotic value 2£ 2 . Thus 



A(K,t) = fU(K)t 



-1/3 



2tt 



-2/3 



which implies, if we identify the scaling parameter with 
the localization length, £(K) ~ £, 

f(x) = 2x 2 . 



Figures fTTTa) and fl2Ta) show the results of the fit- 
ting procedure applied to the numerical data and to the 
experimental data, respectively. In both cases, the proce- 
dure groups all points in a single curve, within the accu- 
racy of the data. The resulting curve clearly displays 
two branches, a diffusive (red) and a localized (blue) 
one, with the critical point being at the tip joining the 
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Figure 12: (Color online) Finite-time scaling applied to the 
experimental results (from 30 to 150 kicks). The scaling pro- 
cedure is the same as in Fig. [TT] (a) The fact that all experi- 
mental points lie on a single curve, with a diffusive (red) and a 
localized (blue) branch, is a proof of the relevance of the one- 
parameter scaling hypothesis, (b) The maximum displayed 
by the scaling parameter £ in the vicinity of K c = 6.4 is a 
clear-cut proof of the Anderson transition. Phase-breaking 
mechanisms (cf. text) smooth the divergence at the critical 
point. When these effects are properly taken into account, 
one obtains a critical exponent v = 1.4 ± 0.3, [the dashed 
line is a fit with Eq. lf48|l ] compatible with the numerical re- 
sults. This plot corresponds to 48 experimental runs. Other 
parameters as in Fig. 1111 



two branches; this is a signature of the Anderson transi- 
tion. It also justifies a posteriori the scaling hypothesis 
Eq. f46|) used for analyzing the data. 

The scaling parameter £(K) is plotted in Figs. fTTTb) 
and fl~2lfb). for numerical and experimental data respec- 
tively. As stated above, this parameter can be identi- 
fied to the localization length in the localized regime. In 
the diffusive regime, it scales as the inverse of the diffu- 
sive constant. Indeed, in the far diffusive regime one has 
(p 2 ) = D(K)t, which implies 



A(K,t) 

m 



D(K)t 



1/3 



so that £(K) = 1/D(K) in the far diffusive regime. 
One notes that £ (K) increases rapidly in the vicinity of 
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the critical value K c , on both sides of the transition. This 
corresponds to a divergence of the localization length and 
to a vanishing of the diffusion constant at criticality (in 
practice smoothed by decoherence, see below). This con- 
stitutes a clear experimental evidence of the Anderson 
phase transition. 

D. Experimental determination of the critical 
exponent 

The behavior of £ (K) gives a fundamental information 
about the transition, namely the value of the localiza- 
tion length critical exponent v. There is a discrepancy in 
the literature between the theoretical predictions v = 1 
[13], v = 1.5 [73] . a nd the result of numerical simulations 
v = 1.57±0.02 [12l | . which stresses even more the impor- 
tance of an experimental determination. In this section, 
we present the first unambiguous experimental determi- 
nation of the critical exponent of the Anderson transition 
in 3 dimensions. 

The finite-time scaling procedure allows us to extract 
from finite-time experimental data the localization length 
I (corresponding to t — > 00), which is the order parameter 
of the Anderson transition. It is given by the scaling 
parameter £(K) and predicted to diverge at criticality 
with the power law 

l~\K-K c \- v . (47) 

We thus expect that the singularity in £(K) can be de- 
scribed by Eq. (j47]l , and to be able to extract the value of 
the critical exponent v. This is of primary importance, as 
there is presently no unambiguous accurate experimen- 
tal determination of v for non-interacting particles, and 
there is a discrepancy in the literature between the theo- 
retical predictions v = 1 [13], v = 1.5 fr3], and the result 
of numerical simulations v = 1.57 ± 0.02 [121 ]. 

When the slope of In A vs lni -1 / 3 is small, as it is near 
the critical point, the scaling procedure tends to round 
off the singularity in £(K). Moreover, decoherence in the 
experiment produces a cut-off the algebraic divergence. 
If the system has a finite phase-coherence time r v , a new 

characteristic length [73] P v = [Dt^] 1 ^ 2 appears in the 
problem, which sets an upper bound for the observable 
localization length I and thus smooths its divergence at 
criticality. In practice, we model such smoothing by in- 
troducing a small cut-off on the divergence of £ (K) , which 
takes into account both the finite-time scaling procedure 
itself and decoherence effects: 

J— = a\K - K c \" + f3 . (48) 

The experimental data have been fitted with this for- 
mula (|48[i [dashed curve in Fig. 112( b)]. which gives K c ~ 
6.4 ± 0.2, and a critical exponent v = 1.4 ± 0.3. In order 
to compare these results to the ideal case of the perfectly 
coherent quasiperiodic kicked rotor, Eq. lj20|) . we also fit- 
ted the curve in Fig. [TTb with Eq. (J48J) ; in this case, the 



cutoff P accounts for limitations of the finite-time scal- 
ing procedure. The model Eq. |48]) fits very well to the 
numerical data [dashed curve in Fig. [TTb ] and gives the 
critical stochasticity K c ~ 6.4±0.1 and the critical expo- 
nent v = 1.6 ±0.2. The good agreement between the nu- 
merical simulations and the experimental results proves 
that spurious effects (such as decoherence) are well under 
control. Moreover, the experimental value we obtained 
v = 1.4±0.3 is compatible with the value found in numer- 
ical simulations of the true random 3D Anderson model 
[HI, [I2]- We emphasize that there are no adjustable pa- 
rameters in our procedure, all parameters are determined 
using the atoms themselves as probes. 



VI. UNIVERSALITY OF THE ANDERSON 
TRANSITION 

At this point, a reasonable question is: Does the 
quasiperiodic kicked rotor exhibits the same critical phe- 
nomena - i.e. belongs to the same (orthogonal) universal- 
ity class [66j] _ as the true 3D-Anderson model. Can this 
simple three-frequency dynamical system exactly mimic 
the critical behavior of 3D disordered electronic conduc- 
tors? In this section, we show that the answer is posi- 
tive: The 3-frequency quasiperiodic kicked rotor and the 
true 3D- Anderson model belong to the same universality 
class. This is a strong claim that relies on a very precise 
determination of the critical exponent v. The accuracy of 
this determination is comparable to that obtained in the 
most sophisticated numerical studies of the 3D Anderson 
model [3,[74j]. Within numerical uncertainties, the criti- 
cal exponent is found to be universal and identical to the 
one found for the 3D- Anderson model [13]. The techni- 
cal details of the calculation have already been reported 
in [73]. We here just discuss the essential ingredients 
proving universality. 



A. Precise estimate of the critical exponent 

Reliably distinguishing the different universality 
classes of the Anderson transition requires a very pre- 
cise determination of the critical exponent; for instance, 
the value v = 1.43 ± 0.04 for the unitary symmetry class 
is close to the one for the orthogonal symmetry class [7(| 
v= 1.57 ±0.02. 

The main uncertainty in our experimental determina- 
tion of the critical exponent is due to statistical errors on 
IIo and to the limited duration of the experiment. How- 
ever, numerical simulations are not limited to 150 kicks 
but can be ran for several thousands of kicks, and sta- 
tistical uncertainties on (p 2 ) can be sharply reduced by 
averaging over initial conditions. The numerical inaccu- 
racy in the finite-time scaling determination of v from 
the numerical data is thus mainly due to the procedure 
failing to reproduce the singular behavior of the scaling 
function at the critical point. 



16 




Figure 13: (Color online) Dynamics of the quasiperiodic 
kicked rotor in the vicinity of the critical regime. The rescaled 
quantity In A(K, t) vs. K is plotted from t = 30 to t = 40000. 
All curves intersect, to a very good approximation, at a single 
point (K c ~ 6.4, lnA c ~ 1.6). This multiple crossing indicates 
the occurrence of the metal-insulator transition. Small devi- 
ations from crossing are due to the existence of an irrelevant 
scaling parameter at finite time and residual correlations in 
the disordered potential (see text). K and e are swept along 
the straight line drawn in Fig. [6] Parameters are £ = 2.85, 

LO'2 = 27T\/5, <^3 = 27T\/l3. 



How can one improve the accuracy on the determina- 
tion of the critical exponent v? This can be achieved 
by fitting directly the raw data lnA(iT, t). The start- 
ing point of our analysis is the behavior of the scaling 
function T = InF in the vicinity of the critical point: 



InA = InF 



(K - K c ) t 



l/3v 



T \(K - K c ) t l ' iv 



(49) 



As lnA(/v, t) is an analytical function for finite t 
(Fig. [TTJ , the scaling function T can be expanded around 



lnA(i) ~ In A c + (K - K c ) t xl2,v T\ + 



(50) 



where In A c = ^[0] and T\ = dJ r (x)/dx| a;= o- 

A remarkable feature of Eq. ([50)) is that when In A is 
plotted against K, the curves for different times t should 
intersect at a common point (A' c ,lnA c ); and this cross- 
ing, indicates the occurrence of the metal-insulator tran- 
sition. This is clearly visible in Fig. [13l Another inter- 
esting feature of Eq. l(50|) is that the critical exponent v 
can be determined from the slope of InA at K c : 



(InA)' (K c ,t) 



<91nA(A", t) 



dK 



oc t 



l/3v 



(51) 



K=K C 



This is the simplest procedure to evaluate the critical ex- 
ponent: (In A)'{K C , t) is evaluated by linear regression of 
In A vs K in a small interval near K c , giving an exponent 
v ~ 1.61 ± 0.10 (see Fig. refng:simulambdaprimevslnt). 
The linear regime has nevertheless very small size: (K — 
Kcjt 1 / 3 " <C 1, and neglecting non-linear corrections lead 




Figure 14: (Color online) Linear regression of (In A)'(i!T c ,t), 
Eq. (HU vs. hit for t = 30 to t = 40000 permits to extract 
the critical exponent v from the slope 1/3^, which is v = 
1.61. It is difficult to assess the uncertainty associated with 
this measurement as it depends crucially on the interval of K 
where the behavior of In A vs K can be assumed to be "linear". 
The parameters are the same as in Fig. 1131 



to systematic errors on the estimation of v. This is why 
the error ±0.1 refers to systematic errors and not to 
the uncertainty on the fitting parameters, which is much 
smaller as easily seen in Fig. [TH 

In practice, there are small systematic deviations from 
Eq. f50|) . Such deviations can have different sources: 

• the presence of an irrelevant scaling variable, that is 
when, in addition to [K — K c )t x l Zv ', (InA) depends 
also on another scaling variable which vanishes in 
the limit t — > oo, but still plays a role at short time; 

• non-linear dependence of the scaling variables in 
the stochasticity parameter K; 

• resonances due to the periods being well approxi- 
mated by a ratio of small integers. 

The latter one is specific to our three-frequency dynam- 
ical system, but the former two also play an important 
role in the standard Anderson model (ill . H3 |. These 
small corrections can be taken into account - following 
the method devised in [l2j for the Anderson model - by 
slightly modifying the basic scaling law, Eq. (J45J) , in two 
ways: introduce a non-linear of the argument of the T 
function with K — K c in Eq. lj49j) on the one hand, and 
allow to subtract irrelevant scaling corrections to (InA) 
on the other hand. To minimize the the effect of res- 
onances, we only retain data for sufficiently long times 
and average over different initial conditions, i.e. different 
quasi-momenta f3 and phases <p% and 953. 

We computed InA for times up to t = 10 6 kicks with 
an accuracy of 0.15%, for which more than 1000 initial 
conditions are required. We analyze data over the full 
range of times t £ [lO 3 , 10 6 ] : The best fit is determined 
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Figure 15: (Color online) lnA s , after subtraction of correc- 
tions due to the irrelevant scaling variable, plotted versus 

In r^/t 1/,3 l and the scaling function deduced from the model 

(black curve). The parameters are that of the set T> (see 
Table [J). The best fit estimates of the critical stochastic- 
ity and the critical exponent are in this case: K c = 8.09 ± 
0.01, In A c = 1.64 ± 0.03 and v = 1.59 ± 0.01. 
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u>3 


K 


E 


A 


2.85 


2tt\/5 


2tt\/I3 


6.24 -> 6.58 


0.413 -> 0.462 


B 


2.85 


2tt\/7 


2tt\/T7 


5.49 -> 5.57 


0.499 -> 0.514 


C 


2.2516 


I/t? 


l/rf 


4.98 -> 5.05 


0.423 -> 0.436 


V 


3.5399 


n/r) 




7.9 -> 8.3 


0.425 0.485 



Table I: The four sets of parameters considered: k, L02 and W3 
control the microscopic details of the disorder, while e drives 
the anisotropy of the hopping amplitudes. 



by minimizing the deviation 



E 

K,t 



\nA(K, t)-F(K,t) 



(52) 



where cr(K, t) is the numerical uncertainty (one standard 
deviation) of the computed quantities hxA(K,t). 

In Fig. HH we plot the scaling function corrected from 
the irrelevant scaling variable, as a function of ^(K)ft 1 ^ 3 . 
All data collapse almost perfectly on the scaling function 
deduced from the model. 

Since the measurement errors in the data introduce 
some uncertainty in the determination of the fitted pa- 
rameters, the confidence intervals for the fitted param- 
eters were estimated using the bootstrap method which 
yields Monte-Carlo estimates of the errors in the fitted 
parameters The fitted parameters presented below 
are given with the corresponding 68.2% confidence inter- 
vals (standard errors). 





K c 


In A c 


V 


y 


A 


6.36 ±0.02 


1.60 ±0.04 


1.58 ±0.01 


0.71 ±0.28 


B 


5.53 ±0.03 


1.08 ±0.09 


1.60 ± 0.03 


0.33 ±0.30 


C 


5.00 ±0.03 


1.19 ±0.15 


1.60 ± 0.02 


0.23 ±0.29 


V 


8.09 ±0.01 


1.64 ±0.03 


1.59 ±0.01 


0.43 ±0.23 



Table II: Best fit estimates of the critical parameters K c 
and In A c , the critical exponent v together with their uncer- 
tainty (one standard deviation), v is expected to be universal 
whereas lnA c and K c do depend on anisotropy [H2| and k, 
LU2 and UJ3. Irrelevant parameters are sensitive to microscopic 
details, therefore y is strictly positive and not universal. 



B. Universality of the critical exponent 

A key property of the Anderson transition is that the 
critical behavior can be described (76l . [zH] in a frame- 
work of universality classes. This means that the critical 
behavior should not be sensitive of the microscopic de- 
tails but should depend only on the underlying symme- 
tries of the system (e.g. time-reversal symmetry). Irrel- 
evant parameters become negligible for sufficiently long 
times/large system size, whereas the relevant parame- 
ter behavior is universal. This brings the universality of 
the critical exponents. When considering a system with 
pseudo-random disorder such as the quasi-periodic kicked 
rotor, one could ask whether the universality is broken 
or not due to correlations in the disorder potential. To 
answer the question, we changed some parameters that 
govern the microscopic details of the disorder potential 
of the quasi-periodic kicked rotor, namely k, UJ2 and W3 
and the path along which we cross the transition. 

The computer time required in those sophisticated nu- 
merical studies is very long. Therefore we chose to re- 
strict ourselves to the detailed study of only four different 
cases, see Table |H 

The estimated critical parameters and their confidence 
intervals are given in Table HTl A typical scaling function 
is drawn in Fig. [TH 

The most important point to be drawn from Table ITT1 
is that the estimates of the exponent v for the four differ- 
ent sets are in almost perfect agreement with each other 
and with the estimate of v based on numerical studies 
of the true random Anderson model v — 1.57 ± 0.02 of 
the orthogonal symmetry class (l2 [. Note also that in 
the case of the quasiperiodic kicked rotor, the critical 
stochasticity K c and lnA c depend on: (i) the anisotropy 
governed by the parameter e and (ii) %, LO2 and 0J3. The 
dependence (i) of the critical disorder and critical In A 
on anisotropy is a typical feature of the Anderson tran- 
sition in anisotropic solids (Hll I52L [53| . The quasiperi- 
odic kicked rotor may indeed be seen to correspond to a 
model of random chains (coupled by terms scaling like e 
in the two transverse directions) considered in [52|], see 
Eq. f33|) . The dependence (ii) follows from the relation 
between the initial "classical" diffusion constant (see sec- 
tion Hll and the parameters %, W2 and w 3 . Such a depen- 
dence was observed both numerically and experimentally 
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for the standard kicked rotor (79l.[80|. and was accounted 
for in terms of correlations between the kicks by Shep- 
elyansky in his early work [79| . 

The Anderson transition with the quasiperiodic kicked 
rotor is a robust feature: we observed that, for certain 
mutually incommensurate triplets (k , uj 2 , u>s) , systematic 
deviations to scaling (such as resonances) can occur for 
intermediate times, but eventually vanish. 

VII. CONCLUSION 

We discussed in detail in the present work the first 
unambiguous evidence of the Anderson transition in 3D 
with atomic matter waves with atomic matter waves 
by realizing experimentally a quasiperiodic kicked rotor. 
This allowed us to put into evidence the existence of the 
transition and to measure its critical exponent thanks 
to a finite-time scaling procedure. Our numerical result 
v = 1.59 ± 0.01 is in perfect agreement with the current 
value for the Anderson model, and is compatible with 
our experimental determination 1.4 ± 0.3. We have also 
shown that the quasiperiodic kicked rotor exhibits the 
same critical phenomena as the truly random Anderson 
model, and therefore that both systems belong to the 
same (orthogonal) universality class. 

These results are particularly relevant since they show 
that it is possible to explore a system like the Ander- 
son model, that played an important hole in many ar- 



eas of physics but resisted thorough experimental inves- 
tigations. One can guess that this kind of analogy will 
be extended to other models in the near future, as evi- 
denced by the work of Wang and Gong [8l| concerning 
the analogy of a quantum kicked rotor and the Harper 
model. This shall open new and exciting tracks in cold- 
atom physics. These analog models can even prove more 
flexible and more powerful than the original ones, as, for 
example, our Anderson-equivalent system can very easily 
be extended to higher dimensions by introducing new in- 
commensurate frequencies. Intermediate situations like 
a 2D kicked rotor with two or three incommensurate 
frequencies might be a convenient solution from the ex- 
perimental point of view. This can hardly be done in 
condensed-matter systems or even in the ultracold atom 
realization of the ID Anderson model [22] • The theo- 
retical study of quantum phase transitions in high di- 
mensions will most probably be boosted as experimental 
results become available. We are presently working in 
this direction: Numerical and experimental determina- 
tions of the critical exponents in four dimensions seems 
feasible. 
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